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Abstract — The problem of definition of effective material 
parameters (permittivity and permeability) for composite layers 
containing only one-two parallel arrays of complex-shaped in- 
clusions is discussed. Such structures are of high importance for 
the design of novel metamaterials, where the realizable layers 
quite often have only one or two layers of particles across 
the sample thickness. Effective parameters which describe the 
averaged induced polarizations are introduced. As an explicit 
example, we develop an analytical model suitable for calculation 
of the effective material parameters e e ff and /u e a for double 
arrays of electrically small electrically polarizable scatterers. 
Electric and magnetic dipole moments induced in the structure 
and the corresponding reflection and transmission coefficients are 
calculated using the local held approach for the normal plane- 
wave incidence, and effective parameters are introduced through 
the averaged fields and polarizations. In the absence of losses both 
material parameters are purely real and satisfy the Kramers- 
Kronig relations and the second law of thermodynamics. We 
compare the analytical results to the simulated and experimental 
results available in the literature. The physical meaning of the 
introduced parameters is discussed in detail. 

Index Terms — Metamaterial, effective medium parameters, 
permittivity, permeability, polarization, local held, causality, pas- 
sivity, reflection, transmission 



I. Introduction 

The problem of extraction of material parameters for com- 
posite slabs implemented using complex-shape inclusions (of- 
ten referred to as metamaterial slabs) has been discussed in 
a large number of recent articles, see e.g. [1]— [10] for some 
example contributions. The physical meaning of many of the 
recent retrieval results is, however, controversial. Typically, the 
5-parameter retrieval procedure based on the Fresnel formulae 
is used to extract the effective material parameters under the 
normal plane-wave incidence. This technique is often known 
to lead to nonphysical material parameters that violate the 
second law of thermodynamics: Either the real part of e e ff 
or ;U e ff obeys antiresonance and therefore the corresponding 
imaginary part has a "wrong" sign, see [3], [7] for some results 
indicating this phenomenon, and [4]-[6] for related criticism. 
It is important to bear in mind that even though the afore- 
mentioned material parameters might satisfy the Kramers- 
Kroning relations (thus being causal), these parameters have 
no physical meaning in electromagnetic sense as physically 
meaningful material parameters must satisfy simultaneously 
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both the causality requirement and the passivity requirement 
[11]. The above discussed parameters (from now on we do not 
use for them the term "material parameter") allow, however, 
to correctly reproduce the scattering results for normal plane 
wave incidence. 

One particular example of complex slabs, referred to in the 
beginning of the paper, is double arrays of small scatterers. 
Such slabs present a very interesting problem as in several 
recent papers the authors state that double grids of electrically 
short wires (or electrically small plates) can be used to produce 
negative index of refraction in the optical [12]— [15] and in the 
microwave regime [16]. In some cases, however, the nonphys- 
ical (antiresonant) behavior is seen in the extracted material 
parameters [14], [16] casting doubts over the meaningfulness 
of assigning negative index of refraction for these structures. 

The goal of this work is to formulate the effective material 
parameters for double arrays of small scatterers through the 
macroscopic polarization and magnetization, and compare 
the results to those obtained using the S-parameter retrieval 
method (the reference results are obtained by different authors 
[14]). The particles in the arrays are modeled as electrically 
small electric dipoles in order to calculate, by using the local 
field approach, currents induced to each of the particles at nor- 
mal plane wave incidence. Knowing the induced currents, the 
reflection and transmission coefficients for the arrays, as well 
as the induced electric and magnetic dipole moment densities 
(macroscopic polarization and magnetization) are calculated. 
Effective material parameters are then defined directly using 
the expressions for the averaged fields and the macroscopic 
polarization and magnetization. Finally, the double array is 
represented as a slab of a homogeneous material characterized 
by effective material parameters. 

Since the layer is electrically thin and spatially dispersive, 
it is obvious that an effective material parameter model 
cannot adequately and fully describe the layer's electromag- 
netic properties, and, moreover, effective parameters can be 
introduced in different ways. It is important to understand 
the physical meaning of different effective medium models 
and their applicability regions. In this paper we compare the 
material parameter extraction based on macroscopic properties 
calculated by means of the local field approach with the 
conventional method based on the 5-parameter inversion. 

The paper is organized in the following way: In Section II 
we use the local field approach to calculate the electric dipole 
moments induced in the particles and define the averaged 
fields needed in the determination of the effective material 
parameters. The derivation of these parameters is presented in 
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Section III. In Section IV some calculated example results are 
presented, and these results are qualitatively compared with the 
numerical and experimental results available in the literature. 
Discussion is conducted on the physics behind the results. The 
work is concluded in Section V. 

II. Local held approach to determine the dipole 

MOMENTS AND THE AVERAGED FIELDS 
A. Electric dipole moments induced in the particles 

Consider a double array of scatterers shown in Fig.Q] At this 
point the shape of the particles can be arbitrary, provided that 
they can be modeled using the electric dipole approximation. 
With this assumption the response of each particle can be 
described in terms of the induced electric dipole moment p 
which is determined by the particle polarizability a and the 
local electric field at the particle position. 

The array is excited by a normally incident plane expressed 
as (see Fig. 02 a )) 



E c 



E e jkoz = E xe jkoz . 
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The common notation for the wavenumber fc 
for the free-space wave impedance rj = J p, / e is used. The 
electric dipole moment induced in a reference particle sitting 
in the first and the second grids reads: 



located in the same grid as the reference particle, whereas 
j3(h) takes into account the influence of the other grid. 

By combining Eqs. (Q~|), © and (O, the following system 
of equations is obtained for the unknown dipole moments: 



a- 1 pi = E + /3(0)pi+/3(/i)p 2 , 
x p 2 = Ece-^ + (3(h)Pi + /3(0)pa. 



(4) 



Solving this system of equations, the electric dipole moments 
induced on the particles in layers 1 and 2 read: 
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where A = [c 



(3(0)] 2 — (3(h) 2 . Approximate analytical 



formulas for the interaction coefficients have been established 
in [17]-[20]: 
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where ~E l °% is the local electric field exciting the reference 
particles. Since both layers are identical in this analysis, the 
notation for the particle polarizability is simplified as ct\ = 
«2 = a- The local electric fields can be expressed as 
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where /3(0) is called the self-interaction coefficient and (3(h) 
is the mutual interaction coefficient. Physically, (3(0) takes 
into account the contribution to the local field of the particles 
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where So is the unit cell area (for the rectangular-cell array 
So = a x b, Fig. OlbX and parameter R is equal to a/1.438 
[18]). The imaginary parts of the interaction constants (3(0) 
and (3(h) are exact, so that the energy conservation law is 
satisfied. The real parts are approximate, and the model is 
rather accurate for ka < 1.5... 2 [18]-[20]. Note that the 
present definition of polarizability a includes also the effects 
of scattering, thus remaining a complex number even in the 
absence of absorption. If there is no absorption in the particles, 
the imaginary part of a reads (e.g., [21]): 
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In the far zone, the reflected field of one grid is a plane 
wave field with the amplitude [21] 



pref _ 



(10) 



The scattered plane-wave field at the reference plane z = 
reads: 



Fig. 1. (a) Geometry of a double array of small scatterers (square patches), 
(b) The unit cell. 
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Finally, the reflection and transmission coefficients can be 
written as: 
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B. Definition of the averaged fields 

The total averaged electric and magnetic fields in the 
vicinity of the array are, by definition, the sum of incident 
field and the averaged scattered field: 
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Under plane-wave excitation each grid creates reflected and 
transmitted plane waves having the amplitudes given by 
Eq. ( [Tol l. Those fields are the scattered fields averaged in 
the transverse plane. In order to average the fields inside the 
slab, the field that is scattered by the first grid and propagates 
towards the second grid, and the field scattered by the second 
grid and propagating towards the first grid are considered, 
see Fig. [2] The averaged electric field inside the slab can be 
expressed as follows: 
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' ' /; ' 1 dz. 
(15) 

After conducting the integration, we find the volume-averaged 
electric field: 



E = e-^/ 2 (E - J g(p 1 +p 2 )) 
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Analogously, the averaged magnetic field can be expressed as 
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Fig. 2. Incident and scattered fields in the problem geometry. 



III. Effective material parameters 

Let us start from the common definition of the constitutive 
parameters of an isotropic homogeneous material: 



D = £ E + P = £o£effE, P = ^, 



B = /i H + M = MoMeffH, M 



V 



(18) 
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Here P and M are the volume-averaged polarization and 
magnetization, p and m are the electric and magnetic dipole 
moments induced in the unit cell, and V is the volume of the 
unit cell. Knowing E, H, P, and M for our particular layer 
we define effective material parameters as: 
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Naturally, these parameters have more limited physical mean- 
ing than the usual parameters of a bulk homogeneous sample, 
but they can be used as a measure for averaged polarizations 
in a thin layer. Let us next find these parameters for our 
particular example grid in terms of its geometry and inclusion 
polarizabilities. 

In order to simplify the forthcoming notations, the sum 
and difference of the electric dipole moments with the use 
of Eqs. (|5j, © can be written in the following manner: 

2e -jk h/2 C0S (k Q h/2) 
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It is important to note that for lossless particles the imaginary 
parts of a'^ 1 and a" -1 can be solved exactly from the power 
balance requirement [19]: 



Im^'- 1 } = ^cos 2 (fc V2), 

DO 

Im{a"- 1 } = !^sin 2 (fc V2). 

<3n 



(25) 
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The use of exact expressions for the imaginary parts is 
very important: Approximate relations would lead to complex 
material parameters of lossless grids, where the imaginary part 
would have no physical meaning. 

The total averaged electric polarization which we need to 
substitute in Eq. < f20b reads 



Pi + P2 

V 



(27) 



In the following, we consider lossless particles with the inverse 
polarizabilities 



= Rc{c/- 1 }+ jlmia'- 1 }, 
= Rc{a"- 1 }+jlm{a"- 1 }. 



The imaginary parts in the above relations satisfy d25l ) and 
Substituting first (ED into d27j, and then d27j into ( 1201 ). 
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the effective permittivity obtained after some mathematical 
manipulations is as follows: 



Seft = 1 + 



koh 



Re{a ^tanj — I -—sin M — 



(30) 

Notice that the permittivity is purely real, as it should be since 
the particles in the arrays are lossless. 

For our system of two identical grids of small particles 
the effective permeability [i e g defined by d20b can be also 
expressed in terms of the induced electric dipole moments. 
Considering one unit cell formed by two particles of length 
I with x-directed dipole moments pi t 2, we can write for the 
currents on the particles (averaged along x) I\ : 2 = jwpxp/l- 
The magnetic moment of this pair of particles (referred to the 
unit cell center) is then 



P2) 
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After inserting first (1221 into d3~TT l. and then (l3TT > into d20l >. 
the effective permeability reads: 
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We can again observe that this quantity is purely real, as it 
should be in the absence of losses. Notice that the procedure 
presented here can easily be generalized to different particle 
geometries. 

As is apparent from the definition, these effective parameters 
measure cell-averaged electric and magnetic polarizations in 
the layer. Can they be used to calculate the reflection and 
transmission coefficients from the layer? This question will 
be considered next using a numerical example. 

IV. Numerical examples 

In this section, the double array of square patches studied in 
[14] is considered as a representative example of comparison 
between the effective parameters introduced here and material 
parameters formally extracted from measured or calculated S- 
parameters. 

The geometry of the double-grid array is the same as 
considered above (Fig. [TJ. The unit cell is characterized by 
the following parameters: Edge lengths of the patches t = I, 
the lattice constants a = b, the metal layer thickness p, and the 
distance between the layers equals h. The dimensions of the 
unit cell considered in [14] are the following: w = I = 300 
nm, t = 20 nm, a = b = 650 nm, p = 20 nm, h — 80 nm. 
Since the double array considered in [14] is targeted to the use 
at THz frequencies, the metal thickness p becomes comparable 
to the dielectric spacer thickness h. The physical thickness of 
the slab is D = h + 2p. However, at THz frequencies one 
must take into account nonzero fields inside metal particles 
(here gold is considered, which at these high frequencies is 
usually characterized by the Drude model (e.g., [11])). For this 
reason an effective thickness of D e ff = 1.5D is chosen for the 
analysis, but we should bear in mind that this is only a rough 
estimation, and an accurate determination of this effective slab 
thickness is very difficult. 



Notice that, as it is stated in [20], the limitation ka < 1.5. ..2 
in the analytical model for the interaction constants (Section II) 
is not critical for these calculations until the lattice resonance 
at ka = 2ir is reached, since the reflection properties near 
resonances are mainly determined by the particle resonances. 
In this case, it means that the model is accurate enough up to 
approximately 460 THz. 

The following rough estimates are used for calculating the 
particle polarizabilities. A known [21] antenna model for short 
strip particles is used by setting the particle width equal to 
the particle length. The polarizability is resonant and can be 
characterized in the microwave regime by an equivalent LC- 
circuit, where C is approximated as the input capacitance of 
a small dipole antenna having length I and radius ro = w/4, 
and L is the inductance of a short strip having length I, width 
w, and thickness t [22]: 

7T(//2) £0 
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1.19 + 0.22—— . (34) 
w + t I J 

(The values of l,w,t are inserted in eq. d34l in millimeters. 
The result is in nano-Henrys.) As mentioned above, in the THz 
range the behavior of metals is different from the microwave 
regime. Following the treatment in [23], the penetration of 
fields inside the particles can be represented by an additional 
parallel capacitance C a dd and a series inductance L a dd [23]: 

EoWt 



L 



add 



SQWtljJt, 



(35) 



(36) 



where ui p is the plasma frequency of considered metal (gold 
in this case) and l e g is the effective particle length. The 
physical clarification of the above parameters is available in 
[23]. Due to the cosine current distribution the effective length 
is approximated as ? e g = 21. Finally, the particle polarizability 
reads: 

= (C + C ddd )l 2 

1-^{L + L M ){C + G M ) 

Fig. [3] shows the magnitudes of the reflection and transmis- 
sion coefficients calculated using Eqs. (fT2l and ([T3l >. Com- 
paring this result with the experimental results presented in 
Fig. 4a of Ref. [14], good agreement is observed. The shape 
of the reflection and transmission functions follows the ones 
presented in [14]. The differences in the amplitude levels and 
transmission frequencies are well expected: The model used 
here does not take losses in the particles into account. 

The effective material parameters calculated using Eqs. ( f30b 
and d32l are depicted in Fig. |U Both material parameters 
clearly behave in a physically sound manner: They are purely 
real, since lossless particles have been considered, and they are 
growing functions outside the resonant region. Contrary, the 
material parameters extracted in [14] for the same structure 
behave in a nonphysical manner: The imaginary parts of 
permittivity and permeability have the opposite signs over 
certain frequency ranges, thus, the passivity requirement is 
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violated [4], [5], [11]. It is evident also from Fig. 4c and 
4d of Ref. [14] that the extracted material parameters do not 
satisfy the Kramers-Kronig relations (e.g., the antiresonance 
in the real part of permittivity cannot correspond to a positive 
imaginary part (with the authors' time dependence assump- 
tion) if inserted into Kramers-Kronig relations), thus, they have 
quite limited physical meaning [11]. Both of these problems 
are avoided with the method described here. On the other 
hand, the material parameters extracted from the reflection 
and transmission coefficients reproduce these properties of 
the layer exactly, while the parameters introduced here do 
not necessarily give so accurate predictions for reflection and 
transmission coefficients. 

This is illustrated next by representing the double layer of 
particles as a slab of a homogeneous material having thickness 
D g and characterized by the material parameters shown in 
Fig. |U The standard transmission-line equations are used to 
calculate the reflection and transmission coefficients for a 
normally incident plane wave. The result is depicted in Fig. [5] 
When comparing Figs. [3] and it can be observed that the 
results do not correspond exactly to each other, even thought 
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Fig. 3. Calculated reflection and transmission coefficients. 
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Fig. 4. Calculated effective permittivity and permeability. 
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Fig. 5. Reflection and transmission coefficients calculated when the double 
array is represented as a slab of homogeneous material characterized by 
e e ff: Uttt depicted in Fig. [4] 

the principal features such as reflection frequencies around 
200 and 300 THz seen in Fig. [3] are reproduced in Fig. [5] 
This is an expected result: There is no such uniform material 
(with physically sound material parameters) that would behave 
exactly as the actual double grid of resonant particles. In 
a truly homogeneous material the unit cell over which the 
averaging is performed should contain a large number of 
"molecules" yet remaining small compared to the wavelength. 
It is clear that for the considered double array this is not the 
case, and also spatial dispersion effects cannot be neglected. 

In particular, very close to the grid resonances the induced 
dipole moments are very large, which corresponds to large 
values of the present effective parameters and to large elec- 
trical thickness of the homogenized slab. These large values 
correctly describe large polarizations in the grids, but they fail 
to correctly predict the reflection and transmission coefficients. 
On the contrary, the effective parameters retrieved from the 
layer's S- parameters, correctly describe these reflection and 
transmission coefficients, but they do not describe polarization 
properties of the actual structure. 

It is also important to note that the standard S'-parameter 
retrieval procedure often leads to nonphysical material pa- 
rameters (see the discussion above). The material parameters 
defined through macroscopic polarization and magnetization 
are always physically sound since they are defined directly 
by the actual dipole moments induced in the microstructure. 
Thus, the effective material parameters assigned in this way 
always make physical sense. 

V. Conclusion 

In this paper an analytical model to assign effective material 
parameters e e ff and /x e g for double arrays of electrically small 
scatterers has been presented. The induced electric dipole mo- 
ments have been calculated using the local field approach for a 
normally incident plane wave, and the corresponding averaged 
fields have been determined. The effective material parameters 
have been defined directly through the macroscopic polariza- 
tion and magnetization. The derived expressions have been 
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validated by comparing the calculations to the numerical and 
experimental results available in the literature. The main con- 
clusion from this study is that simple effective medium models 
of electrically thin layers with complex internal structures 
are always limited in their applicability and in their physical 
meaning. Some models are best suitable to describe reflection 
and transmission coefficients at normal incidence (like the 
parameters conventionally retrieved from S-parameters), other 
models describe well the averaged induced polarizations in 
the structure and allow one to make conclusions about, for 
instance, negative permeability property (like the parameters 
introduced in this paper). In applications, it is important to 
understand what model is used and what properties of the layer 
this particular model actually describes. Finally, it is necessary 
to stress that the present study has been restricted to a particu- 
lar special case of a dual layer of planar electrically polarizable 
particles. This approach needs appropriate modifications and 
extensions if, for instance, inclusions are also magnetically 
polarizable. For slabs containing more than two layers, the 
direct extension of this model corresponds to averaging over a 
unit cell containing three or more particles, which apparently 
is not useful for electrically thick slabs. 
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